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ABSTRACT 

We propose a model to explain a puzzling 3:2 frequency ratio of high frequency 
quasi-periodic oscillations (HFQPOs) in black hole (BH) X-ray binaries, GRO J1655- 
40, GRS 1915-1-105 and XTE J1550-564. In our model a non-axisymmetric magnetic 
coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc co- 
exists with the Blandford-Znajek (BZ) process. The upper frequency is fitted by a 
rotating hotspot near the inner edge of the disc, which is produced by the energy 
transferred from the BH to the disc, and the lower frequency is fitted by another ro- 
tating hotspot somewhere away from the inner edge of the disc, which arises from the 
screw instability of the magnetic field on the disc. It turns out that the 3:2 frequency 
ratio of HFQPOs in these X-ray binaries could be well fitted to the observational data 
with a much narrower range of the BH spin. In addition, the spectral properties of 
HFQPOs are discussed. The correlation of HFQPOs with jets from microquasars is 
contained naturally in our model. 
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1 INTRODUCTION 

Quasi-periodic oscillations in X-ray binaries have become a 
very active research field since the launch of the Rossi X- 
Ray Timing Explorer {RXTE; Bradt, Rothschild & Swank 
1993). A key feature in these sources is that some of high 
frequency quasi-periodic oscillations (HFQPOs) appear in 
pairs. Five black hole (BH) X-ray binaries exhibit transient 
HFQPOs, of which three sources have pairs occurring in 
GRO J1655-40 (450, 300Hz; Strohmayer 2001a, hereafter 
sola; Remillard et al 1999), GRS 1915-^105 (168, 113Hz; 
McChntock & Remillard 2003, hereafter MR03), and XTE 
J1550-564 (276, 184Hz; Miller 2001; Remillard et al 2002) 
with a puzzling 3:2 ratio of the upper frequency to the lower 
frequency. Very recently, the 3:2 frequency ratio (henceforth 
3:2 ratio) was discovered in the X-ray outburst of H1743- 
322 (240, 160Hz; Roman et al. 2004; Remillard et al. 2004). 
These discoveries give the exciting prospect of determining 
BH mass and spin, as well as testing general relativity in the 
strong-field regime. 

A number of different mechanisms have been proposed 
to explain the origin of HFQPO pairs in BH X-ray bina- 



ries. Strohmayer (SOla; 2001b) investigated combinations of 
the azimuthal and radial coordinate frequencies in general 
relativity to explain the HFQPO pairs in GRO J 1655-40 
and GRS 1915-1-105. Wagoner et al. (2001) regarded the 
HFQPO pairs as fundamental g-mode and c-mode discoseis- 
mic oscillations in a relativistic accretion disc. Abramowicz 
& Kluzniak (2001) explained the pairs in GRO J 1655-40 as 
a resonance between orbital and epicyclic motion of accret- 
ing matter. It seems that more than one physical model is 
required to fit all of the HFQPO observations. 

As is well known, the Blandford-Znajek (BZ) process 
is an effective mechanism for powering jets from AGNs and 
quasars, and it is also applicable to jet production in mi- 
croquasars (Blandford & Znajek 1977, Mirabel & Rodrigues 
1999). Recently, the magnetic coupling (MC) of a rotating 
BH with its surrounding disc has been investigated by some 
authors (Blandford 1999; Li 2000a, 2002a; Wang et al. 2002, 
hereafter W02), which can be regarded as one of the variants 
of the BZ process. The MC process can be used to explain a 
very steep emissivity in the inner region of the disc, which is 
found by the recent XMM-Newton observation of the nearby 
bright Seyfert 1 galaxy MCG-6-30-15 (Wilms et al. 2001; Li 
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2002b; Wang et al. 2003a, hereafter W03a). In addition, we 
explained the HFQPOs in X-ray binaries based on a model 
of a rotating hotspot due to the MC of a rotating BH with 
the inner region of the disc (Wang et al. 2003b, hereafter 
W03b). However, the 3:2 ratio of HFQPOs has not been 
discussed by virtue of the MC process. 

Very recently, we discussed the condition for the coex- 
istence of the BZ and MC processes, and found that this 
coexistence always accompanies the screw instability of the 
magnetic field in BH magnetosphere, provided that the BH 
spin and the power-law index for the variation of the mag- 
netic field on the disc are greater than some critical values 
(Wang et al. 2003c, 2004, hereafter W03c and W04, respec- 
tively). 

In this paper we propose a model to explain the 3:2 
ratio of HFQPOs in BH X-ray binaries by virtue of non- 
axisymmetric MC of a rotating BH with its surrounding 
accretion disc. The upper frequency arises from a rotating 
hotspot near the inner edge of the disc, and the lower fre- 
quency is produced by the screw instability somewhere away 
from the inner edge. It turns out that the 3:2 ratio of HFQ- 
POs in GRO J1655-40, GRS 1915+105 and XTE J1550-564 
is well fitted by energy transferred in the MC process, where 
six parameters are used to describe the BH mass, spin and 
distribution of a non-axisymmetric magnetic field. The 3:2 
ratio obtained in our model provides a much narrower range 
of the BH spin compared with the other models. Further- 
more, the jets from these BH binaries are explained naturally 
by the BZ process coexisting with the MC process. 

This paper is organized as follows. In section 2 we give 
a brief description of our model. In section 3 we fit the up- 
per and lower frequencies of HFQPOs by the inner and outer 
hotspots co-rotating with the disc. The discussion is focused 
at the issues of the outer hotspot. A scenario of how the 
presence of the screw instability leads to the outer hotspot 
is given. In addition, the time-scale of the screw instabil- 
ity is estimated by using an equivalent R-L circuit, and the 
spectral properties of HFQPOs are discussed by assuming 
the existence of corona above the disc. In section 4 we dis- 
cuss some issues of astrophysical meanings related to the 3:2 
ratio. Throughout this paper the geometric units G = c = 1 
are used. 



2 DESCRIPTION OF OUR MODEL 

We intend to fit the 3:2 ratio of HFQPOs by virtue of an ro- 
tating BH surrounded by a magnetized accretion disc based 
on the following assumptions. 

Q The configuration of the poloidal magnetic field is 
shown in Fig. 1. The radius Vms indicates the inner edge of 
the disc, being the innermost stable circular orbit (ISCO) 
(Novikov & Thorne 1973, hereafter NT73), and is the 
radius where the screw instability of the magnetic field might 
occur, which is related to the angle 6s by the same mapping 
relation given in W04. As shown in Fig. 1, 9s is the angular 
boundary between the open and closed field lines on the 
horizon, and is the lower boundary angle for the closed 
field lines connecting the BH with the disc. Throughout this 
paper 9l ~ 0.457r is assumed in calculations. 

^ The poloidal magnetic field on the BH horizon is 



Astrophysical Load 




Figure 1. Poloidal magnetic field connecting a rotating BH with 
a remote astrophysical load and the surrounding disc. 



non-axisymmetric as assumed in W03b, varying with the 
azimuthal angle (ji as follows (see Fig. 2 in W03b), 



B^(0) = J((B^)^)/(0), /(0) 



1 + S, < 0< A0, 
1, ^ s; 2tt, 



where y {{5^)'^) is root- mean-square of the poloidal mag- 
netic field over the angular coordinate from 9 — to 6l- The 
parameter S is used to describe the strength of the bulging 
magnetic field in the azimuthal region < < A(j}. 

Q The magnetosphere is assumed to be force-free out- 
side the BH and the disc, and the closed magnetic field lines 
are frozen in the disc. The disc is thin and perfectly con- 
ducting, lies in the equatorial plane of the BH with the inner 
boundary being at ISCO. 

@ The poloidal magnetic field varies as a power law 
with the radial coordinate of the disc as follows. 



(2) 



where represents the poloidal magnetic field on the disc, 
and the parameter n is the power-law index, and the dimen- 
sionless radial coordinate is defined as ^ = r/rms- The radius 
rms is related to the BH mass M and the dimensionless spin 
parameter a* = a/M as given in NT73. 

@ The magnetic flux connecting the BH with its sur- 
rounding disc takes precedence over that connecting the BH 
with the remote load. As argued in W03c and W04, as- 
sumption 5 is crucial for the coexistence of the BZ and MC 
processes. 

The magnetosphere anchored in a Kerr BH and its sur- 
rounding disc is described in Boyer-Lindquist coordinates, 
in which the following parameters are involved (MacDonald 
and Thorne 1982, hereafter MT82). 

= (r^ + a'Y - A sin^ 9, p'^ = + cos^ 9., 
A = -I- - 2Mr, vj = (E/p) sin 9, (3) 



In W03b the BZ and MC powers in non-axisymmetic 
case are related to those in axisymmetic case by 

(4) 



(1) 
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MC 



(5) 



where Pg^ and Pmc ^"^^ the powers in non-axisymmetic 
case, and Pg^ and Pmc ^^e the powers in axisymmetic case, 
respectively. The parameter A is expressed in terms of 5 and 
e = A(j)/2n by 

X^[{l + S)e + {l-e)f = {l + Sef . (6) 

Similarly, the BZ and MC torques in non-axisymmetic 
case are related to those in axisymmetic case by 



T, 



BZ 

NA 
MC 



= AT, 



flZ) 
MCi 



(7) 

(8) 



„4 /„ 9 / fc (1 — fc) sin 
P^z/Po^2at I ^ ' 



Pmc/Po = 2al 



2-(l-g) sin-^e» 

Pjl- P) sin^ Ode 
2 - (1 - g) sin^ 6 ' 



o<e <es, 



(9) 



9s<e<9L. (10) 



Tiz/To 



4a, (1 

^MC /To = 
4a, (1 



e (1-13) sin^ Ode 
2-{l-q) Bin-' ' 



o<e <es, 



(11) 



2-(l-,)sin^9 ' es<e<ei,. 

In equations (|^ — 1121 1. k and /3 are the ratios of the angular 
velocities of the magnetic field lines to the angular velocity 
of the BH horizon in the BZ and MC processes, respectively. 
The quantity q = -^/l — a» is a function of the BH spin, and 
Po and To are defined as 



Po^{{B^^f)M' 
To = UBlf)M'': 



• Bird^BH X 6.59 x 10^** erg • s"\ 
Blm^H X 3.26 x lO'^^g ■ crn^ ■ s" 



where B4, is the strength of the poloidal magnetic field on 
the horizon in units of IQ'^gauss, and m^^ is defined as 
ruBH = M/Mq. 

So far six parameters are involved in our model, in which 
rrigjj and a, are used for describing the mass and spin of the 
Kerr BH, and the three parameters, B4, 5 and e are used 
for describing the non-axisymmetric magnetic field on the 
horizon, the power-law index n is for the variation of the 
magnetic field with the radial coordinate of the disc. 



3 INNER AND OUTER HOTSPOTS FOR 3:2 
FREQUENCY RATIO 

As argued in W03b, the upper frequency of HFQPOs is 
modulated by an inner hotspot, which corresponds to the 
maximum of function Fqpq as follows. 



Fqpo = r^F/rl,,Fo = ^""{Fda + Fmc)/Fo, 



(14) 



where Fda is the radiation flux due to disc accretion, and 
Fmc is radiation flux arising from the energy transferred 
from the rotating BH into the inner disc in the MC process. 
We find that the radial coordinate of the inner hotspot, ^max, 



only depends on the parameters a,, and n, i.e., it is inde- 
pendent of the parameters m^^, B4, 5 and e. Considering 
that the hotspot is frozen at the disc, we have the upper 
frequency i'q'pq by substituting ^max into the Keplerian 
angular velocity as follows, 

(15) 



ti-i/2 3 I \ 



where T^i and T^i^ are the torques in non-axisymmetic 
case, and Tj^^ s-^d 2^mc the torques in axisymmetic case, 
respectively. The powers Ppz a-nd Pmc^ and the torques 
Tg^ and T^jq are given in W02 as follows. 



where uo = (msi/)"^ x3.23x 10'*/fz. It turns out that Vq'po 
depends on the three parameters, m^jj, a* and n. 

The lower frequency of HFQPOs is fitted by an outer 
hotspot rotating with the Keplerian angular velocity ex- 
pressed by equation 115II . and a scenario for the production 
of the outer hotspot is given as follows. 

It is well known that the magnetic field configurations 
with both poloidal and toroidal components can be screw in- 
stable (Kadomtsev 1966; Bateman 1978). According to the 
Kruskal-Shafranov criterion (Kadomtsev 1966), the screw 
instability will occur, if the toroidal magnetic field becomes 
so strong that the magnetic field line turns around itself 
about once. Recently, some authors discussed the screw in- 
stability in BH magnetosphere (Gruzinov 1999a; Li 2000b; 
Tomimatsu et el. 2001). In W04 we suggested that the cri- 
terion for the screw instability in the MC process could be 
expressed by 



Fsc 



(a, ; 5, n) — L/{2tvujd) < 0, 



(16) 



where L is the poloidal length of the closed field line con- 
necting the BH and the disc, zu^ is the cylindrical radius on 
the disc, and Fscrew (a»; ^) is a function of a,, ^ and n. 

Setting equality in equation 1161 . we obtain the criti- 
cal radial coordinate = rs/vms for the screw instability, 
which depends on the parameters a, and n. 

It is believed that a disc is probably surrounded by 
a high-temperature corona analogous to the solar corona 
(Liang & Price 1977; Haardt 1991; Zhang et al 2000). Very 
recently, some authors argued that the coronal heating in 
some stars including the Sun is probably related to dissipa- 
'■^"^■^tion of currents, and very strong X-ray emissions arise from 
variation of magnetic fields (Galsgaard & Parnell 2004; Pe- 
ter et al. 2004). 

Analogously, if the corona exists above the disc in our 
model, we expect that the corona above might be heated 
by the induced current due to the screw instability of the 
non-axisymmetric magnetic field. Therefore a very strong 
X-ray emission would be produced above the rotating outer 
hotspot, and the lower frequency Vqpq is obtained directly 
by substituting into equation 1151 . 

The time-scale of the screw instability can be estimated 
by using an equivalent circuit as shown in Fig. 2, in which the 
segments LM and KN represent the two adjacent magnetic 
surfaces connecting the BH horizon and the outer hotspot as 
shown in Fig. 2a. Considering the existence of the toroidal 
magnetic field threading each loop, we introduce an inductor 
into the circuit, which is represented by the symbol AL in 
Fig. 2b. 



A£H = (A**'/27r)nj/, Aei 



(A*^/27r) Q,D, 



(17) 



where A^''' is the flux of the poloidal magnetic field sand- 
wiched by the magnetic surfaces LM and KN . The resis- 
tance at the BH horizon is given by 

a; 2phMh 



ARh = an 



2TVCUI- 



ZUh 



(18) 
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where an — 4,11 = 377ohm is the surface resistivity of the 
BH horizon (MT82), and A9h is the angle on the horizon 
spanned by the two surfaces. The quantities and ro^ are 
the Kerr metric parameters at the horizon. The inductance 
AL in the R-L circuit is defined by 



AL = A*^/r, 



(19) 



where and A^*"^ are the poloidal current fiowing in the 
circuit and the flux of the toroidal magnetic field threading 
the circuit, respectively. The flux A^*-^ can be integrated 
over the region KLMN as follows, 



A*' 



B ^grrgoedrdO, 



(20) 



loop 



where the toroidal magnetic field measured by "zero- 
angular-momentum observers" (ZAMOs) is given as 



(21) 



where a is the lapse function defined in equation Q. 

Unfortunately, the geometric shapes of the magnetic 
surfaces LM and KN are unknown. As an approximate 
estimation we assume that the surfaces LM and KN are 
formed by rotating the two adjacent circles around the sym- 
metric axis, and A^^ can be calculated by integrating over 
the region KLMN as shown in Fig. 3. Finally, we calcu- 
late the quantity r = AL/ARh by incorporating equations 
1181 — 1211 1. which is independent of 7^. 

As a simple analysis we can divide one event of the screw 
instability into two stages. In stage 1 the instability starts, 
and the magnetic energy is released rapidly, the poloidal 
current and thus the toroidal magnetic field decrease to zero 
in very short time. In stage 2 the toroidal magnetic field 
recovers gradually with the increasing poloidal current in 
the R-L circuit, until the next event of the instability occurs 
according to the criterion 11611 . It seems reasonable that the 
duration in stage 1 is much less than that in stage 2, and 
the time-scale of the screw instability depends mainly on the 
duration of stage 2, which is regarded as the relax time of 
the poloidal current in the R-L circuit. The current 
is governed by the following equation, 



dl^ 

AL— + ARhI^ = Ash + Ael- 
at 

Setting the initial condition, — 



(22) 



0, we have the solution, 
I''{t)=ll,^^,,{l-e-''-), (23) 
is the steady current in the R-L circuit, and it 



where P, 

stei 

reads 



I^,,,,y = {AeH + A£L)/AR„. (24) 

From equation 12311 we obtain that the poloidal current 
attains 99.3% of I^t^^^y in the relax time treiax = 5r, im- 
plying the recovery of the toroidal magnetic field. Thus the 
time-scale of the screw instability can be estimated as 

iScrew ^ trelax — 5t. (25) 

Summarizing the above results, we have the upper and 
lower frequencies of HFQPOs, and the corresponding val- 
ues of ^max, Cs and treiax as showu lu Table 1. It is easy to 
check from Table 1 that the time-scale of the screw insta- 
bility is generally greater than the corresponding periods of 
HFQPOs, i.e., 




AL 



Ae„ 



AR„ 



(b) 



Figure 2. (a) Two adjacent magnetic surfaces connecting the BH 
horizon and the outer hotspot, (b) one loop of equivalent circuit 
for screw instability. 



i^Screw ^ f^relax ^ / QPO / Q 



lower 
PO ■ 



(26) 



As pointed out in MR03, the spectral properties of ac- 
creting black holes in states which show HFQPOs are usually 
dominated by a steep power-law (SPL) component. How- 
ever, the origin of X-ray power-law remains controversial. 
Most of models for the SPL state invoke inverse Compton 
scattering as the operant radiation mechanism, and the ra- 
diation spectrum formed through Comptonization of low- 
frequency photons in a hot thermal plasma cloud may be 
described by a power law (Pozdnyakov, Sobol & Syunyaev 
1983). It is noted that the origin of the Comptonizing elec- 
trons involves magnetic instabilities in the accretion disc 
(Poutanen & Fabian 1999), with which the screw instability 
in our model is consistent. 

Suppose that the spectral properties of the inner and 
outer hotspots arise from blackbody spectra, the effective 
radiation temperature (Ths)^jj is expressed by 



Teff — To [/da + /mc]''^'^ , 

To = {Fo/af^^ ^ 4.8 x lO'sy" (K) , 



(27) 
(28) 
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Table 1. The 3:2 ratio of HFQPOs produced by the inner and outer hotspots with S = 0.5 and e = 0.2. Henceforth BHB-I, H and HI 
represent GRO J1655-40, GRS 1915+105 and XTE J1550-564, respectively. 



BHB 


a. 


n 




^relax (sGc) 


Inner Hotspot 




Outer Hotspot 




^max 


jpupper / r>l/2 

(kev) 


upper 

'^QPO 
(Hz) 


?sc 


jT^lower / d1/2 
^HS / ^4 

(kev) 


..lower 
'^QPO 
(Hz) 


I 


0.777 


6.22 


6.8 


7.61 X 10-=* 


1.277 


7.23 X 10-3 


450 


2.015 


4.96 X 10-3 


300 




0.730 


6.27 


5.8 


6.57 X 10--^ 


1.281 


6.72 X lO--^' 


450 


2.097 


4.46 X 10-3 


300 


II 


0.773 


6.23 


18 


2.02 X 10-2 


1.277 


7.19 X 10-3 


168 


2.021 


4.92 X 10-3 


113 


0.603 


6.17 


10 


1.11 X 10"^ 


1.297 


5.55 X 10--' 


168 


2.428 


3.27 X 10-3 


113 


III 


0.788 


6.19 


11.5 


1.28 X 10-2 


1.276 


7.37 X 10-3 


276 


2.000 


5.09 X 10-3 


184 


0.696 


6.29 


8.5 


9.72 X 10--^ 


1.284 


6.38 X 10-3 


276 


2.159 


4.14 X 10-3 


184 



Note: The value ranges of the BH mass corresponding to GRO J1655-40, GRS 1915+105 and XTE J1550-564 are adopted from 
Greene et al. (2001), MR03 and Orosz et al. (2002), respectively. 



4 DISCUSSION 



Inspecting Table 1, we find that the 3:2 ratios of HFQPOs 
in GRO J1655-40, GRS 1915+105 and XTE J1550-564 are 
weU fitted by adjusting the two parameters, a* and n, for the 
given value ranges of m^^ . In this section we discuss some 
issues of astrophysical meanings related to the 3:2 ratio of 
HFQPOs in the BH binaries. 











M 


\ . 1 N 


n B i 


/ Accretion/Disc 



Black HoleX 



Figure 3. Two adjacent magnetic surfaces produced by two ro- 
tating circles, where the arrows represent the current flowing in 
the R-L circuit. 



where a is the Stefan-Boltzmann constant. The energy of 
the inner and outer hotspots can be estimated by 

Ehs = kBT.ff = Eo [foA + juct'^ , 



Eo = kB {Fo/af^ 



4.14 X 10"^By^ 



(kev) 



(29) 
(30) 



where ks is the Boltzmann constant. Substituting ^max and 
into equation 12911 . we have the energy of the inner and 
outer hotspots £;J|'='' and £;jf^='' as shown in Table 1. 

Inspecting equation 1291 and the expressions for fo.A 
and /mc given in W03b, we find that both E^g" and 
^iomer (ijgpgjjfj jjq^ Quly ou the parameters, a* and n, but 
also on the parameters for the non-axisymmetric magnetic 
field, B\^^ , S and e. It is worthy to point out that both 
E-JT and E'^T 
X-ray binaries. 



are independent of the BH mass of the 



4.1 3:2 ratio and magnetic field on the BH 
horizon 

It is pointed out that the twin HFQPOs are not always 
detected together in GRO J1655-40, i.e., 300 Hz QPO is 
observed sometimes without 450 Hz detection, while 450 
Hz QPO is observed sometimes without 300 Hz detection 
(sola). In our model the upper and lower frequencies of 
HFQPOs are produced by the inner and outer hotspots, 
which are located at different sites of the disc. In addition, 
our model is proposed based on the non-axisymmetric mag- 
netic field on the horizon, and two different physical mecha- 
nisms are involved, i.e., the MC process for the inner hotspot 
and the screw instability for the outer hotspot. Since the 
non-axisymmetric magnetic field cannot be stationary on 
the horizon, we can understand that the twin HFQPOs do 
not have to occur simultaneously in some cases. 

In Table 1, both E^"^'"' and Sjf^"' reach the energy 
level of emitting X-ray, provided that the root-mean square 
of the magnetic field on the horizon is strong enough to be 
~ 10^gau.s.s {B4 ~ 10^). Considering that a system cannot 
radiate a given fiux at less than the blackbody tempera- 
ture (Frank et al. 1992), the strength of magnetic field of 
~ lO^gauss might be regarded as a upper limit for energy 
level of emitting X-ray in BH binaries. 

It is pointed out in SOla that the energy spectra of the 
300 Hz QPO appear to be significantly different from those 
of the 450 Hz QPO in GRO J1655-40: The former has a 
typical amplitude in the 2-12 keV band, while the latter is 
detected in the hard band, 13-27 keV. It is also shown in 
Fig. 4.16 of MR03 that the power densities corresponding to 
the upper frequency are greater than those corresponding to 
the lower frequency of HFQPOs in some BH X-ray binaries, 
such as GRO J1655-40 (450, 300Hz) and XTE J1550-564 
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(276, 184Hz). These observations are consistent with the 
results obtained in our model, i.e., E'^g'^^ is greater than 
j^i-o-wer £^j. gach BH binary as shown in Table 1. 

4.2 BH spin and mass constrained by 3:2 ratio 

It is widely believed that HFQPOs in BH binaries might 
be a unique timing signature constraining the BH mass and 
spin via a model rooted in general relativity. Thus a precise 
measurement of the 3:2 ratio of HFQPOs becomes an ap- 
proach to estimate the BH mass and spin. Combining the 
fitting results in Table 1 with the observations of the three 
BH binaries, we discuss the issues related to the BH mass 
and spin as follows. 

The BH spin in GRO J1655-40 has been estimated by 
some authors. Interpreting the X-ray spectra, Zhang, Cui & 
Chen (1997) suggested that the most likely value is 0.7 < 
a, < 1, while Sobczak et al. (1999) give an upper limit a* < 
0.7. Assuming that a bright spot appears near the radius 
of the maximal proper radiation flux from a disc around a 
rotating BH, Gruzinov (1999b) inferred an upper limit a, < 
0.6 for the BH in GRO J1655-40. Abramowicz & Kluzniak 
(2001) suggested that the 3:2 ratio in GRO J1655-40 arises 
from a resonance between orbital and epicyclic motion of 
accreting matter, and they estimated that the value range 
of the BH spin is 0.2 < a, < 0.65. Compared with the 
above estimation ranges, we provide a much narrower range 
for the BH spin in GRO J1655-40, 0.730 < a, < 0.777, as 
shown in Table 1, where the BH spins in GRS 19154-105 and 
XTE J1550-564 are also constrained in rather narrow ranges, 
0.603 < a* < 0.773 and 0.696 < a. < 0.788, respectively. 

4.3 3:2 ratio and jets from microquasars 

The BH X-ray binaries, GRO J1655-40, GRS 1915-^105 and 
XTE J1550-564, are also regarded as microquasars, from 
which relativistic jets are observed (Mirabel & Rodrigues 
1998, 1999 and references therein; see also Table 4.2 of 
MR03). This correlation can be explained very naturally by 
virtue of the magnetic field configuration as shown in Fig. 1. 
The jet power is assumed to be provided by the BZ mech- 
anism and expressed by equations (|1J and Substituting 
a,, n, nigg and into these equations, we have the jet 
powers as shown in Table 2. In calculations the parameter 
k = 0.5 is taken for the optimal BZ power (MT82). 

Mirabel & Rodriguez (1999) pointed out that GRS 
1915+105 may have a short-term jet power of ~ l{!p^erg/ s 
in its very high state, which is a large fraction of the observed 
accretion power. From Table 2 we have the angular bound- 
ary and the corresponding BZ power for GRS 1915-1-105 as 
follows. 

Ps/ = 4.87 X 10^^ (ery/s) , < 6) < 50.5", (31) 

P^/ = 4.68 X 10^** (ery/s) , < 6(5 < 42.4°. (32) 

Equations 13111 and 13211 imply that a large fraction of 
the rotating energy is extracted from the BH contained in 
GRS 1915-1-105 for the jet power, and Pbz can reaches ~ 
10^® erg/s, provided that the magnetic field on the horizon 
is strong enough to attain ~ l(f gauss {B4 ~ 10^). This 
strength of the magnetic field is the same as that required 
by the inner and outer hotspots for emitting X-ray. By using 
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n 
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a. 



(b) 

Figure 4. The BZ power (a) varying with n with the given 

values of a, , (b) varying with a» with the given values of n. 

equations Q S'Hd @ we find that the BZ power increases 
monotonically with the power-law index n for the given BH 
spin a,, while it varies non-monotonically with a, for the 
given n as shown in Fig. 4. 

4.4 3:2 ratio and electric current on disc 

As shown in Table 1, a very large power-law index n varying 
from 6.17 to 6.29 is required by the 3:2 ratio of HFQPOs in 
the BH binaries. It implies that the poloidal magnetic field 
increases inward in a very steep way governed by equation 

It is worthy to point out that this result is consistent 
with that in W03a, where the power-law index n varying 
from 6 to 8 is required to explain a very steep emissivity 
(4-5) observed in Seyfert 1 galaxy MCG-6-30-15. 

The poloidal magnetic field arises most probably from 
the toroidal electric current on the disc, and the current 
density j^p is related to by Ampere's law as follows. 

.^^J_dB^^ _J_dBl ^ _n(B^^_(„+,) ^^^^ 
47r dr 47rrms 47rr-ms 

Inspecting equation 1331 . we obtain two features of the 
toroidal electric current on the disc. 
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Table 2. The jet powers accompanying 3:2 ratio of HFQPOs with 5 = 0.5 and s = 0.2. 



BHB 




n 




ds 


P^i/Blierg/s) 


I 


0.777 


6.22 


6.8 


50.7" 


7.09 X 10^8 


0.730 


6.27 


5.8 


48.9" 


3.98 X 10^" 


II 


0.773 


6.23 


18 


50.5" 


4.87 X 10^9 


0.603 


6.17 


10 


42.4" 


4.68 X 10^" 


III 


0.788 


6.19 


11.5 


50.9" 


2.14 X 10^9 


0.696 


6.29 


8.5 


47.7" 


7.01 X 10^« 



(i) The minus sign in equation 1331 shows that the cur- 
rent flows in the opposite direction to the disc matter; (ii) 
the power-law index n -I- 1 in equation 1331 implies that the 
profile of the current density in the central disc is steeper 
than that of the magnetic field. We give a primary explana- 
tion for these features as follows. 

Probably the direction of the current arises from the 
effect of the closed magnetic field lines on the charged par- 
ticles of the disc matter. In our model the angular velocity 
of the BH is required to be greater than that of the disc to 
produce the inner and outer hotspots on the disc. Therefore 
the charged particles of the disc matter will be dragged for- 
ward by the closed field lines. Considering that the mass of 
the positive charged particle, such as proton, is much greater 
than that of electron, we think that the bulk fiow of elec- 
trons might be a little ahead of that of protons. Probably the 
minute difference in the two kinds of the bulk flows results 
in a current in the opposite direction to the disc. 

On the other hand, the flux of angular momentum 
transferred from the rotating BH to the disc is given by 



rrNA 



4:Tvr 



(34) 



Substituting equations (|HJ and 1121 1 into equation 13411 . we 
have 



Incorporating equation |(2J with equation H35|l . we have 
Hmc oc C • 



(35) 



(36) 



Considering that the flux of angular momentum, H^jq, 
transferred to the disc will block the accreting particles, and 
the blocking effect would be much more signiflcant for elec- 
trons than protons due to the mass differences. Thus more 
electrons would stay in the region very close to ISCO. That is 
why the current density consist of the bulk flow of electrons 
and it varies with ^ in a power law expressed by equation 

As a summary we intend to give a comment for the six 
parameters involved in our model. The BH mass m^^ has 
its range given by observations, and B4 ~ 10"" is required by 
the hotspots for emitting X-ray. The parameters e and 5 are 
used to describe the non-axisymmetric magnetic fleld on the 
BH horizon. In fact, neither the position nor the energy of 
the inner and outer hotspots is sensitive to the values of e 
and S, and we take 5 = 0.5 and e = 0.2 in calculations. Thus 
the rest two parameters, a, and n are adjustable in fitting 



the 3:2 ratio, which play very important roles in fitting the 
3:2 ratio of HFQPOs. 

In this paper a non-axisymmetric magnetic field on the 
BH horizon is assumed in an ad hoc way. Unfortunately, we 
cannot give a good explanation for its existence at present. 
Probably the non-axisymmetric magnetic field on the BH 
horizon arises from non-axisymmetric accretion disc, since 
the magnetic field is brought and held by the surrounding 
magnetized disc. We hope to approach this difficult task in 
future. 
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